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1. Introduction 

Heavy quarkonia (charmonium, bottomonium, ...) provide an ideal set of observ- 
(~| , ables to probe properties of low-energy QCD in a controlled way. The reason is the 

^ ' following. Heavy quarkonia are non-relativistic bound states and, therefore, char- 

acterized by a set of energy scales hierarchically ordered: m, mv, mv^, ... where m 
is the heavy-quark mass and v <^ 1 the relative heavy-quark velocity. For heavy 
' quarkonia, m is much larger than the scale of non-perturbative physics, Aqcd, and, 

therefore, degrees of freedom associated with that scale can be treated perturba- 
tively and calculations done order by order in as. The non-relativistic hierarchy 
of scales also survives below Aqcd- Therefore, for any heavy quarkonium state the 
low-energy dynamics is organized in matrix elements ordered in powers of v (and, in 
general, Aqcd/™)- To any given order in as and v, only a finite number of Feynman 
diagrams and matrix elements respectively have to be calculated. 

The way to implement rigorously these expansions in QCD is provided by the 
non-relativistic effective field theories (EFTs) of QCD. The first has been Non- 
Relativistic QCD, NRQCEr-i^. It is obtained from QCD by integrating out degrees 
of freedom of energy m. NRQCD still contains the lower energy scales as dynamical 
degrees of freedom. In the last few years, the problem of integrating out the remain- 
ing dynamical scales has been addressed by several groups and has now reached a 
solid level of understanding (lists of references may be found i iP). The ultimate 
EFT obtained by subsequent matchings from QCD, where only the lightest de grees 
of freedom of energy mv^ are left dynamical, is potential NRQCD, pNRQClJ^'^l 
This EFT is close to a quantum-mechanical description of the bound system and. 
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therefore, as si mple . It has been systematically explore d in the dy namical regime 
mv'^ > Aqcd in^'^ and in the regime mu^ ^ Aqcd iiJMIlalailiil, An alternative 

approach to pNRQCD has been suggested i This approach has been devel- 
oped in the dynamical situation mv'^ ^ Aqcd, but not been extended, so far, to 
mv'^ ^Aqcd where most of the heavy quarkonium states are believed to lie. 

In this letter, I will review the theory status of heavy quarkonium inclusive and 
electromagnetic decays into light particles in the framework of QCD non-relativistic 
EFTs. The main mechanism of heavy quarkonium decay into light particles is quark- 
anti-quark annihilation. Since this happens at a scale 2m, which is perturbative, the 
heavy quarks annihilate into the minimal number of gluons allowed by symmetry. 
Experimentally this fact is reflected by the narrow width of the heavy quarkonia 
below the open flavour threshold. In an EFT language, once the scale m has been 
integrated out, the information on decays is carried by contact terms (four-fermion 
operators) whose matching coefficients develop an imaginary part. The low-energy 
dynamics is in the matrix elements of the four-fermion operators evaluated on the 
heavy quarkonium states. If one assumes that only heavy quarkonium states with 
quark-anti-quark in a singlet configuration can exist, then only singlet four-fermion 
operators contribute and the matrix elements reduce to heavy quarkonium wave 
functions (or derivatives of them) calculated in the origin. This assumption is known 
as the "singlet model". Explicit calculations show that at higher order the singlet 
matching coefficients develop infrared divergences (for P waves this happens at next- 
to-leading ordeili^ compare with the expressions in [Appendix A| ) . In the singlet 
model, these do not cancel in the expression of the decay widths. It has been the 
first success of NRQCD to show that, due to the non-Abelian nature of QCD, 
the Fock space of a heavy quarkonium state may contain a small component of 
quark-anti-quark in an octet configuration bound with some gluonic degrees of 
freedom (the component is small because operators coupling transverse gluons with 
quarks are suppressed by powers of v), due to this component, matrix elements 
of octet four-fermion operators contribute and, finally, exactly these contributions 
absorbe the infrared divergences o f the singlet matching coefficients in the decay 
widths, giving rise to finite result J^^^l. NRQCD is now the standard framework 
to study heavy quarkonium decays. From the theoretical side in recent years the 
main effort has gone into two obvious directions: (1) improving the knowledge of 
the perturbative series of the matching coefficients either by fixed order calculations 
or by resumming large contributions (renormalons or large logs); (2) improving the 
knowledge of the NRQCD matrix elements either by direct evaluation, which may 
be obtained by fitting the experimental data, by lattice calculations, and by models, 
or by exploiting the hierarchy of scales still entangled in NRQCD and constructing 
EFTs of lower energy. I already mentioned that pNRQCD is the ultimate of these 
EFTs. In such context a new factorization can be achieved that allows to reduce, 
under some dynamical circumstances, the number of non-perturbative parameters 
in the expression of the decay widths. 
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Experimentally several facilities have operated and are operating in an energy 
range relevant for heavy quarko niu m. I refer to the pages of the Quarkonium Work- 
ing Group for a broad overvie Here, I would like to mention only some of the 
data produced in the last few years relevant for heavy quarkonium inclusive and elec- 
tromagnetic decays. They come from the E835 experiment at FNAL, where heavy 
quarkonium is produced from pp annihilation (E835 operated in the charmonium 
energy region) and from the B-factories (BABAR at SLAG and BELLE at KEKB, 
which operate at the T(4S') resonance), BES at BEPG (operating in the charmo- 
nium energy region) and GLEG at GESR (GLEOIII took data in the bottomonium 
energy region, GLEO-c is taking data in the charmonium one) where heavy quarko- 
nium is produced from e^e~ collisions. In the bottomonium system, GLEO has pro- 
vided the first experimental values (still affected by large uncertainties) for the ratios 
of the inclusive decay widths of 2Pj bottomonium states, extracted from the data 
on two-photons transitions from T(3S') decays In the charmonium system, new 
determi nat ions of the rjc resonance parameters came from the experiments E SSd^, 
BABAEfni 

(where t he ri r. is produced from tw o-ph oton interactions), BELLB I I 
{B — > rjc K) and BEy^i*' {J/i^ ^ Vcl)-^ refer td^for comparison and discussion of 
the data. E835 also provides T{r]c — > 77^^- The resonance parameters of the •0(25') 
have been newly extracted from high statistics e+e~ cross-section data at BEs!22. 
The same collaboration also reports a new determination of r('0(2S') e+e^). The 
branching ratio for ip{2S) — > e~^e~ has been measured by the E835 experimenu2Sl 
The branching ratios of •0 (25' ) — > e"'"e~ and ■0(25') — > have been measured 

by the BABAR experimenl!^. BELLE reported the first observation of the r^(25) 
and the measurement of its resonance parameters (from B — > r]c{'2S) K irM^ and 
e+e" ^ J/ipr]c{2S) followed by BABAR^^and GLEcEi (in both cases from 

two-photon interactions). For what concerns the L = 1 charmonium states, the 
Xco resonance parameters have been newly meas ured at E83^^. The sa me c ollab- 
oration provides a determination of T{xcQ 7 7 j'^^^'^' ^ and T{xc2 t -tI^ . New 
values for T{xcO ~^ u) and r(xc2 ^77) have been provided by GLEO^^and for 
r(Xc2 ^77) by BELLEP2in both cases from two-photon production processes. 

The letter is organized as follows. In Sec. |21 I review the NRQGD factorization 
formulas for heavy quarkonium decay widths. In Sec. 01 1 briefly discuss the pertur- 
bative series of the matching coefficients. [Appendix Al contains a complete list of all 
the imaginary parts of the matching coefhcients of the dimension 6 and 8 operators 
(hadronic and electromagnetic) at their present accuracy. In Sec. 0] I discuss the 
NRQGD matrix elements and in Sec. 14.11 the pNRQCD factorization. Some final 
remarks are given in Sec. El 



2. NRQCD 

The NRQGD factorization formulas are obtained by separating contributions com- 
ing from degrees of freedom of energy m from those coming from degrees of freedom 
of lower energy. In the case of heavy quarkonium decay widths, they have been rig- 
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orously proveo^^ High-energy contributions are encoded into the imaginary parts 
of the four-fermion matching coefficients, f , gi,s,.ee,-yi(^^^^ Lj) and are ordered in 
powers of as ■ Low-energy contributions are encoded into the matrix elements of the 
four-fermion operators on the heavy quarkonium states \H) ((. . .)h = {H\ . . . \H)). 
These are, in general, non-perturbative objects, which can scale as powers of Aqcd, 
mv, TOw^, ... (i.e. of the low-energy dynamical scales of NRQCD). Therefore, ma- 
trix elements of higher dimensionality are suppressed by powers of v or Aqcd/'ti- 
Including up to the NRQCD four-fermion operators of dimension 8, the NRQCD fac- 
torization formulas for inclusive decay widths of heavy quarkonia into light hadrons 
(LH) reacPSl 

+lmfsCSi) {OsCSi))vQ(nS)+lmfsCSo) {OsC So))vQ(nS) 



T{PQ{nS) -^LH) = ^ |^Im/i('5o) (Oi (^5o))p«(„s) 

V{xQ{nJS) -^LH)^^ flm ^^i?'^' Pj)) x.inJS) 

+lmh{^'+'Ss) (08('5o))x«(nJS)^ ■ (3) 

At the same order the electromagnetic decay widths are given by: 

2 



T{VQ{nS) ^ e+e-) = ^ (^Im/ee^Si) {O^mC Si))vQ(nS) 



T{PQinS) ^ = A ( ini/^^(i5o) {OemC So)) P^^nS) 



^Such a proof is still lacking for the NRQCD factorization of heavy quarkonium production cross 
sections. 
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(7'em(^5'o))pq(„s) 



(5) 



r(xQ(nJl)^77) = 2Im^^(3p; 



{OEM{^Pj))xQ{nJl) 



J = 0,2. 



(6) 



The symbols Vg and Pq indicate respectively the vector and pseudoscalar S'-wave 
heavy quarkonium and the symbol XQ the generic P-wave quarkonium (the states 
XQ(nlO) and xq(»^>/1) are usually called hQ{{n — 1)P) and XQj{{''^ ~ 1)^); respec- 
tively). 

The operators O, 'Pi^g„^M{^^^^ Lj) are the dimension 6 and 8 four-fermion oper- 
ators of the NRQCD Lagrangian. They are classified in dependence of their trans- 
formation properties under colour as singlets (1) and octets (8) and under spin (5), 
orbital (L) and total angular momentum (J). The operators with the subscript EM 
are the singlet operators projected on the QCD vacuum. The explicit expressions 
of the operators can be found irJ^ (or listed in [Appendix A| oi'^). 



3. The perturbative expansion 

The imaginary parts of the four-fermion matching coefficients have been calculated 
over the past twenty years by different authors and to different levels of precision. 
Since the results are scattered over a large number of papers, some of them being 
difficult to collect, some having been corrected in subsequent publications and some 
still being in disagreement with each other, I have listed all the imaginary parts of 
the matching coefficients of the dimension 6 and 8 operators (hadronic and electro- 
magnetic) at the present accuracy in [Appendix A| The tree-level matching of the 
dimension 10 S'-wave operators can be found in^. The tree-level matching of the 
dimension 9 electromagnetic P-wave operators can be found i 

The convergence of the perturbative series of the four-fermion matching coeffi- 
cients is often bad. Let us consider, for instance, the following matching coefficients 
(n/ = 3, = 2mf^. 

13.7^), 

as 



ImfsCSo) 
Im/sfSi) 



Im/i('^o) 
Im/i(3p2) 



Apart from the case of Im/i 



1 + 11.1 ■ 



1 + 10.3- 



13.6-0.44 log — 
2 m 



1 - 0.73+ 1.67 log 



TT 



2m/ TT 

^P^), the series in of the other coefficients does not 
show convergence. This behaviour cannot be adjusted by a suitable choice of the 
factorization scale /i, which enters only in Im/i('^Po^2)- A solution may be provided 
by the resummation of the large contributions in the perturbative series coming from 
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bubble-chain diagrams. This analysis has been successfully carried out for S'-wave 
annihilation decayJ^. A treatment that includes P-wave decays is still missing. 

4. The relativistic expansion 

The NRQCD matrix elements may be fitted on the experimental decay datcP^EHIor 
calculated on the latticd^. The matrix elements of singlet operators can be linked 
at leading order to the Schrodinger wave f unc tions in the origirPl and therefore may 
be evaluated by means of potential models^. In general, however, NRQCD matrix 
elements, in particular of higher dimensionality, are poorly known or completely 
unknown. 

It has been discussed iiPancP, that higher-order operators, not included in 
the formulas of Sec. [3 even if parametrically suppressed, may turn out to give 
sizeable contributions to the decay widths. This may be the case, in particular, 
for charmonium, where v'^ ~ 0.3, so that relativistic corrections are large, and for 
P-wave decays where the above formulas provide, indeed, only the leading-order 
contribution in the velocity expansion. In fact it was pointed out in"^^ '^^ that if no 
special cancellations among the matrix elements occur, then the order relativistic 
corrections to the electromagnetic decays Xco — > 77 and Xc2 77 rnay be as large 
as the leading terms. 

In'^' it was also noted that the numerical relevance of higher-order matrix ele- 
ments may be enhanced by the multiplying matching coefficients. This is, indeed, 
the case for the decay width of S'-wave vector states, where the matching coefficients 
multiplying the octet matrix elements (with the only exception of Im/8("^Pi)) are 
enhanced by as with respect to the coefficient Im/i('^S'i) of the leading singlet 
matrix element (see Eq. and [Appendix A| ). 

In the bottomonium system, 14 S- and P-wave states lie below the open flavour 
threshold (T{nS) and rif){nS) with n = 1,2,3; hb{nP) and Xbj{nP) with n = 1,2 
and J = 0, 1, 2) and in the charmonium system 8 {^{nS) and ric{nS) with n = 1, 2; 
hc{lP) and Xcj(l-P) with J = 0,1,2). For these states Eqs. ([ij-® describe the 
decay widths into light hadrons and into photons or e+e" in terms of 46 NRQCD 
matrix elements (40 for the S-wave decays and 6 for the P-wave decays). More 
matrix elements are needed if, as discussed above, higher-order operators have to 
be included. 

4.1. pNRQCD 

The number of non-perturbative parameters may be reduced by integrating out 
from NRQCD degrees of freedom with energy lower than to, since each degree of 
freedom that is integrated out leads to a new factorization. Eventually, one ends up 
with pNRQCD, where only degrees of freedom of energy tow^ are left dynamical. In 
the context of pNRQCD, the NRQCD four-fcrmion matrix elements can be written 
either as convolutions of Coulomb amplitudes with non-local correlators (in the 
dynamical situation mv^ > Aqcd) or products of wave- functions in the origin by 
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non-local correlators (in the dynamical situation mv^ <C Aqcd)- The first si tuation 
is believed to be the relevant one at least for the bottomonium ground statJ^^^EHl 
In the limiting case mu^ » Aqcd, the correlators reduce to local condensates 
and explicit formulas have been worked out for the electromagnetic decay of the 
T(IS') i The last situation is expected to be the relevant one for most of the 
existing heavy quarkonia (wit h the po ssible exception of the bottomonium ground 
state) and has been studied iJHEEl However, a general consensus on the above 
assignations of heavy quarkonium states to dynamical regions has not been reached 
yet. As an example, I mention that iiJ^ it is suggested that also some of the 
higher bottomoniu m stat es may be Coulombic bound states while in practically 
all potential modelJ^SEHl 

the bottomonium ground state is described by means of 

confining potentials. 

In the situation mv^ <C Aqcd and under the condition that: (a) all higher gluonic 
excitations between the two heavy quarks develop a mass gap of order Aqcd, (b) 
threshold effects are small, and (c) contributions coming from virtual pairs of quark- 
antiquark with three momentum of order y^toAqcd are subleading, the NRQCD 
octet matrix elements relevant for Eqs. Ql-lUll can be written at leading order in 
the V and Aqcd I'm expansion a a^'^t 



{Of,{^Si))vQ(nS) - {Os,C So)) PQ(nS) 

27r \ 3m 



— ^ — ^ ' (") 



In^QW (08(^'5'l))pQ(ng) 



(08(3Pj))v^(„5) 



3 

_^ |fii°J(0)P ( {Ca/2~Cf)cIB, \ 
2^ j' 

(08(lPl)}p^(„S) 

^(2J + l)C.^%l^(- (^-/%^-)^^ ), (9) 

(O8(^^0))x«(„J5) = ^ '^"l^^")'' g3, (10) 

' 3 7rm^ 

where cp stands for the chromomagnetic matching coefficient, which is known at 
next-to-leading ordeiPSI Therefore, at the considered order, the octet matrix ele- 

''Concerning the perturbative calculation of the electromagnetic decay width of the T(1S) a renor- 
malization group improved ex press ion can be found in^' and the wave function in the origin at 
next-to-next-to-leading order i 

'^Conditions (a) and (b) select the simplest version of pNRQCD with only one degree of freedom: 
the heavy quarkonium singlet field. Condition (a) is supported bv l attice data on the excitation 
spectrum of the gluon field around a static quark-antiquark paiil^3 Condition (b) may be prob- 
lematic for the ^{2S), whose mass is very close to the DD nrodiiction threshold. Condition (c) is 
more technical and affects the matching to NRQCD. I refer tcjl^lfor a discussion of its validity. 
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ments factorize into the product of the zeroth-order radial part of the heavy quarko- 
nium wave function, i?^*^' , which may be calculated from the real part of the pN- 
RQCD Hamiltoniary, and some chromoelectric and chromomagnetic correlators: 

1 r°° 

fn = ^y^ d""(5E(0-5E(0)), (11) 
I3n = ^ rdte{gB{t)-gBiO)), (12) 

Jo 

1 f°° Z"*^ f*^ f 
AN '^^^ '^^^ dt3(t2-i3)'j ({gE(ti)-,gE(i2)}{ffE(t3)-,5E(0)}), 

- A (Tr(gE(ti) • 5E(i2)) Tr(.gE(<3) • 5E(0)))e|, (13) 



where 

(.gE(ti) • g^{t2) g^{h) ■ 5E(0))e = (5E(ti) • gE(t2) 5E(i3) • <?E(0)) 

- ^(<7E(ti) •gE(t2))(.9E(i3) •5E(0)). (14) 

These correlators are universal in the sense that they do not depend on the heavy 
quarkonium state and , hence, may be calculated once for ev er, e ither by means of 
lattice simulations'^ or specific models of the QCD vacuurrp2l or extracted from 

o 

some set of experimental data . 

At leading order in the v and Aqcd I'm expansion the singlet matrix elements 
can be expressed in terms of the wave functions in the origin onljIS 

{Oi{^Si))vQ{nS) = {OlC Sq)) p^i^nS) = (OeM (^S*! )) („S) 



|i?i°J(0)| 



2 



= (Oem(^S'o))pq(„s) = Ca — , (15) 

zvr 

(Ol('^+'P/))x<,(nJ5) = (OEM('^+^Pj))x«(n,75) = ^^|i?iT(0)r (16) 

At leading order the matrix elements of the Vi operators involve also the correlator 



{'Pl{^Si))vQ{nS) - {'PlCSo))pQ{nS) - CPemC Si))vQ{nS) 



(17) 



where E^°^ ~ M- 2m is the leading-order binding energy. Equation p7l) reduces to 
the formula obtained i if the heavy quarkonium state satisfies also the condition 
mv > Aqcd- 

The leading corrections to the above formulas come from quark-antiquark pairs 
of three momentum of order -^toAqcd- The existence of this degree of freedom 
in the heavy quarkonium system has been pointed out i where the leading 
correction to Eq. (|15|) has been calculated. 
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The pNRQCD factorization formulas reduce the number of non-perturbative pa- 
rameters needed in order to describe heavy quarkonium decay widths'^. In particu- 
lar, they have been used to calculate bottomonium matrix elements from charmo- 
nium data. This is useful since, at the moment, bottomonium data are less abundant 
than charmonium o nes. I n this way P-wave bottomoniu m in clusive decay widths 
have been calculatecPS3 before the first data by CLEO-IlffSl were made available. 
One has to stress, however, that the theoretical uncertainties associated to P-wave 
heavy quarkonium decays are rather large, due to the large corrections either in 
the perturbative series (as discussed in Sec. or in the relativistic expansion (as 
discussed in Sec. . For the inclusive decay width of P-wave heavy quarkonium 
neither the resummation of large perturbative corrections, nor the computation of 
operators appearing at next-to-leading order in the v and Aqcd/"J expansion has 
been done yet. 

5. Conclusions 

In this letter, I have reviewed some general aspects of the theory of inclusive heavy 
quarkonium decays. The standard framework is provided by NRQCD and more 
generally by non-relativistic EFTs of QCD. These have put the study of heavy 
quarkonium observables on the solid ground of QCD. Models and phenomenolog- 
ical approaches have not necessarily become obsolete: they may provide estimates 
of the non-perturbative paramaters that appear in the EFTs. In particular, poten- 
tial models may still be useful to estimate the heavy quarkonium wave functions. 
However, also potentials are parameters of the EFT and have a precise expression 
in terms of the original degrees of freedom (gluons and quarks) of QCD. Lattice 
gauge theories provide the most natural and well founded tool to calculate non- 
perturbative quantities. In fact several lattice determinations of matrix elements 
entering in the heavy quarkonium decay width expression at the level of NRQCD, 
as well as of correlators and Wilson loops entering at the level of pNRQCD already 
exist. 

Experimentally, heavy quarkonium decay data have been produced in large 
amount in the last years and have improved the accuracy of several of the mea- 
sured widths and branching ratios. They call for comparable precise theoretical 
determinations. The relevance is twofold. On one hand we may extract from heavy 
quarkonium data several of the non-perturbative parameters that characterize the 
low-energy dynamics of QCD. This is possible, because we have simple and exact 
expressions that factorize the non-perturbative physics. As an example, I mention 
that the correlators entering in the expression of the decay widths in pNRQCD give 
information on the masses of the heavy quarkonium exotic hybrid states, describe 
the behaviour of the QCD static potential at intermediate distances and contribute 
to the heavy quarkonium levels. On the other hand we may use heavy quarkonium 
data to extract some of the fundamental parameters of the Standard Model (e.g. the 
heavy-quark masses and as). In the case of as, this is not yet possible from heavy 
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quarkonium decay data with an accuracy comparable with other determinations 
due to the difficulties discussed in this review'^. This is one of the many challenges 
in the present and future of heavy quarkonium physics. 



Appendix A. Imaginary parts of the dimension 6 and dimension 8 
four-fermion matching coefficients 

The imaginary parts of the four-fermion matching coefficients are known to different 
levels of precision. In the following, 1 will indicate their most updated values. The 
symbols stand for: Ca = = 3, Cp = (N^ - l)/(2iVc) = 4/3, Bp = (Ca/2 - 
Cf){CI/2-2) = 5/12, Tf = 1/2, /3o = 11 iVc/3-4TFn//3, eg the electrical charge 
of the quark Q (cb = —1/3, Cc = 2/3, ...), a the electromagnetic coupling constant, 
tts the strong coupling constant in the MS scheme and n/ the number of active 
light flavours (typically 3 and 4 for charmonium and bottomonium respectively). 
The scale /i is the factorization scale and the scale the renormalization scale. 
In a physical quantity, like the decay width, the dependence will be canceled by 
low-energy matrix elements and the /i/j dependence by higher-order terms in the 
perturbative expansion. The strong coupling constant Q!s(/ifl) has to be understood 
as running with Uf flavours :° 

'Ca 
2 

/199 13 



Cf 



Cf TTasipLRY 



1 



Cf 



18 24 



Ca 



16 



n/Tp 4- /3o log 



2m 



(A.l) 



9 



Im/i (351)13 1[^^-9)Cf{CI-A){^-Cf] a^{mf 



X <! IH 

IT 



-9.46(2) Cf + 4.13(17) Ca - 1.161(2) + ^/Jq log ^ 

2 TO 



13 „ cks 
— Of — 

4 TT 



(A.2) 



<^The value quoted by the PDcjSl seems to underestimate some of the uncertaintiej^l 
''Note that from 



= f 1 + - V — + ^("')) ' 

V TT 2 ^Jfl J 

'(m) it follows that 

as {^^R) = as ' {^lR) [1 - - —Tplog h 0{a^ ) 

\ i -K m J 



(nf), . (rif+l), \ . r 11 

and Og (m) = (m) it follows that 
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ImA(PO'-y-'^^(^— -C^ja.(^-- + -7r - log — j 



(A.3) 



X <! IH 

IT 



/427 TT^ 



27 



Ini/i(3Po) 



EOI 



3Cf 



X <i iH 

TT 



6 
2 

2 



2m y 2to 
Cf I irasinB,)'^ 



27 



29 

loff 

6 ^ 2m 



■ /3o H 



2m 



(A.4) 



(A.5) 



4 , M 

log 

3 ^ 2to 



(A.6) 



f a. 













— ^^ + -log2 Ca 



29 M 

log 

15 2to 



/3o loj 



Mi? 
2to 



X + 



/2239 337 , 5 

-4Cf + + - 

V 216 384 3 



^- + ^log2 Ca 



- log + /3a log — 
15 2™ / 2m 



(A.7) 



(A.8) 



Inigi(^5o) 



ACf [Ca 



Cp I TTOg , 



(A.9) 
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Imgii^Si) ^- — Cf(C^-4) 



Cf] a! 



X <^ 1 + 



^ 7r2\ ^ /122 17 2 , ^ 



--rvT^+/3olog — 



luifsi-'Si) n/ 



f a„ 













133 2, ^ 7r2\^ 10 ^ 

h - log 2 \Ca nfTp 

18 3^ 4/ 9^ 



73 67 2 \ 5 



36 



. /3olog„ 
n / 2m 



Im/si Pi) = ^™s: 



Im/8(3Po)'^ 3BF7ra,(/i«)' 
X n+ — 



B- 










29 


- log 


2m 



463 35 , ^ 17 2 1 ^ 

^ log 2 TT^ Ca 

81 27 ^ 216 ' 



- /3o log 



2m 



56 



1369 23 



108 18 



Ca + 



Ini/8('P2) ^BF7ras(Mi?)' 
5 



X 1 



4955 43 2 7 



431 72 



^- + ^log2 



29 
15 
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1 



(A.17) 



3 
79 
18 



1 -4Cf 



TT^ - -TT^ log — + 27r2 log 2 



39 



89 2 2, M 5 2, „ 151 13^ 22 , ^ir, ^ 



_ + _iog_)r^„, 



4 , 44 , 

-9- +yi^- 



(A.18) 



Im/,,,(35i) 13 4^e|a3 {l - 9.46(2)C^^} 



(A.19) 



Im/,,(3Po) 1^ i^e%a^ ^ 1 + ^ ( - o + " ) 



TT V 3 4 



(A.20) 



Im/,,(3P2)»1 J.e^a^{l-4C.^}, 



(A.21) 



151 



4 4 2 



(A.22) 



Inigee(^5'i) 



m 



2 2 

-Treoa , 



(A.23) 



(A.24) 



The number over the equal sign indicates the reference/references where the most 
updated value of the matching coefficient can be found. Some comments are in or- 
der. The order as corrections to Im/i('^S'i) and Im /^-^^(^S'l), given in Eqs. I|A.2|I 
and l|A.19p respectively, are known only numerically and, therefore, affected by a 
numerical error. The last line of Eq. I|A.2|) . proportional to a^, comes from the 
annihilation of the quark-antiquark pair into a virtual photon, which then decays 
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into light hadrons. For the order corrections to Im/i(^Po) and Im/i('^P2) there 
are at the moment two (numericaUy shghtly) different determinations in the Uter- 
ature. Since the question of which one, if any, is correct has not been settled yet 
I have reported both determinations (Eqs. (|A.4|I - I|A.5|I and Eqs^ (IA.7|l - l|A.8(l 'l. The 
expression of Im/gf^Sj) given in Eq. (|A.12|I and taken fronJ^ is different from 
that one reported irP^I (there n/ — *■ n-//3). According to one of the authors this is 
the correct onJ^. The electromagnetic matching coefficients (jA.17|) - (|A.24|) refer to 
annihilation processes where the final states consist of two photons (Eqs. (|A.17|I . 
(IA.20p . HA. 211) and HA.22P '). three photons (Eq. (|A.19p ) and 2 massless fermions 
(Eqs. (|A.18I) . HA.23|I and (|A.24p '). Note that the matching coefficient Im/ee(^S'i), 
having been calculated at order a^, is the most accurately known. 

The running equations for the imaginary parts of the matching coefficients of 
the four-fcrmion NRQCD operators of dimension 6 and 8 have been obtained i 
and can be read there in Appendix C. 
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